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Abstract 

The flux jump dynamics in the flux flow regime of type II superconductors is investigated, analyticaUy. It is 
found that under some conditions flux jump avalanche may occur in a superconductor sample, which taJses 
into account an incrtial properties of the vortex matter. 
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Introduction 

As we know, the flux jumps results in a large-scale flux 
avalanches in a superconductor and their origin are related to 
the magnetothermal instabilities [1-5] . Thermomagnetic instabil- 
ity or flux jump phenomena have been observed in conventional 
hard superconductor, as well as in high-tcmperature supercon- 
ductors, recently [1-6]. The spatial and temporal development 
of this instability depends on the sample geometry, temperature, 
external magnetic field, its rate of change and orientation, initial 
and boundary conditions, etc. The critical state instabilities re- 
sult in flux redistribution towards the equilibrium state and are 
accompanied by a signiflcant heat release, which often leads to 
the superconductor-to-normal-transition. Recently, Chabanenko 
et al. [6] have reported an interesting phenomenon in their ex- 
periments - convergent oscillations of the magnetic flux arising 
from flux jump avalanches [6-11]. The authors argued that the 
observed oscillations due to flux avalanches can be interpreted as 
a result of the existence of a definite value of the effective vortex 
mass [12-21]. Thus, it is necessary to take into account collective 
modes, i.e., the incrtial properties of the vortices in studying the 
dynamics of the fiux avalanches. In the present work, we study 
the dynamics of the magnetic flux avalanche, which takes into 
account an inertial properties of the vortex matter. 

1. Formulation 

Bean [1] has proposed the critical state model which is success- 
fully used to describe magnetic properties of typo II supercon- 
ductors. According to this model, the distribution of the mag- 
netic flux density B and the transport current density j inside a 
superconductor is given by a solution of the equation 



^ rotB = jioj. (1) 

(!^^ When the penetrated magnetic flux changes with time, an elec- 
ts trie fleld E{r, t) is generated inside the sample according to Fara- 
^ day's law 
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In the flux flow regime the electric field E(r, t) induced by the 
moving vortices is related with the local current density j{r,t) 
by the nonlinear Ohm's law 

E = vB. (3) 

To obtain qualitative results, we use a classical equation of mo- 
tion of a vortex, which it can derived by integrating over the 



microscopic degrees of freedom, leaving only macroscopic forces 
[21]. Thus, the equation of the vortex motion under the action 
of the Lorentz, pinning, and viscosity forces can be presented as 



dV 

m— + riV + Fl + Fp 



0. 



(4) 



Here m is the vortex mass per unit length, Fl = -J'I'o is the 

c 

Lorentz force, F„ = -jc^o, V = — % — — is the flux flow viscos- 

c C^pn 

ity coefficient, $o ~ vr/ic/2e is the magnetic flux quantum, Hc2 
is the upper critical fleld of superconductor, pn is the normal 
state resistivity, jc is the critical current density [4]. For sim- 
plicity we have neglected the Magnus force, assuming that it is 
much smaller then the viscous force (for example, for Nb sec, 
[6]). In the absence of external currents and flelds, the Lorentz 
force results from currents associated with vortices trapped in 
the sample. 

2. Basic equation 

In combining the relation (3) with Maxwell's equation (2), 
we obtain a nonlinear diffusion equation for the magnetic flux 
induction B{r, t) in the following form 

1 

-I- nl/ = — 

dt 
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(5) 
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The temperature distribution in superconductor is governed by 
the heat conduction diffusion equation 

dT 



v{T) 



dt 



V[n{T)VT]+jE, 



(7) 



Here = J'(r) and k = i^{T) are the specific heat and ther- 
mal conductivity, respectively. The above equations should be 
supplemented by a current-voltage characteristics of supercon- 
ductors, which has the form 

] = j,{T,Bj). 

In order to obtain analytical results of equations (5)-(7), we 
suggest that jc is independent on magnetic field induction B 
and use the Bean critical state model jc = jc{Be,T), i.e., 
jc{T) = jo — a{T — To) [4]; where Be is the external applied mag- 
netic field induction, a = jo/ (Tc — To), To and Tc are the equilib- 
rium and critical temperatures of the sample, respectively, jo is 
the equilibrium current density. For the sake of simplifying of the 
calculations, we perform our calculations on the assumption of 
negligibly small heating and assume that the temperature profile 
is a constant within the across sample and thermal conductivity 



1 



K and heat capacity v are independent on the temperature profile 
[5]. 

We study the evolution of the thermal and electromagnetic 
penetration process in a simple geometry - superconducting semi- 
infinitive sample 2: > 0. We assume that the external magnetic 
field induction is parallel to the z-axis and the magnetic field 
sweep rate is constant. When the magnetic field with the flux 
density Be is applied in the direction of the z-axis, the transport 
current j{x, t) and the electric field E{x, t) are induced inside the 
slab along the y-axis. For this geometry the spatial and temporal 
evolution of thermal and magnetic field perturbations 



Be + b{x,t), 
V = VQ + v{x,t) 



B 



(8) 



are described by the following system of differential equations [8, 
11] 
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(10) 

(11) 



where To{x), Be{x) and Vo{x) are solutions to the unperturbed 
equations, which can be obtained within a quasi-stationary ap- 
proximation. Here we have introduced the following dimension- 
less variables 



B 
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and parameters 
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Here L = cBe/fiojc is the magnetic field penetration depth, 
which is determined from the following equation 



B{x, t) = Be + ^.ojc{x - L), 
with the appropriate boundary conditions 

dB(0,t) = Be, B{L,t) = Q. 

3. Dispersion relation 



(12) 
(13) 



Assuming that the small thermal and magnetic perturbations has 
the form Q{x,t),b{x,t),v{x,t) ~ exp[7f], where 7 is the eigen- 
value of the problem to be determined, we obtain from equations 
(9)-(ll) the following dispersion relations to determine an eigen- 
values of the problem 



(7 + /3) 
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2/3]^ + [(m + 1)7' + 



[(7 + ,fl)M 

- 1]7 - (/i - l)/^lb 



The instability of the flux front is defined by the positive value of 
the rate increase Re 7>0. An analysis of the dispersion relation 
shows that, the grows rate is positive Re 7>0, if > /ic = 2 and 
any small perturbations will grow with time. For the case when 
fi < fic, the growth rate is negative and the small perturbations 
will decay. At the critical value of = fic, the growth rate is 
zero 7=0. For the specific case, where n = 1 the growth rate 
is determined by a stability parameter /3. Thus, the stability 
criterion can be written as 

/3> 1. 

For the case, where thermal effects is negligible (/3 — 1) we may 
obtain the following dispersion relation 



d^b 
dx^ 



■m£ + (7-1)(/^ + 1)6^ 



0. 



(14) 



Seeking for b ~ exp(ifca;) in dispersion relation, the growth rate 
7 dependence can be obtained as a functions of wave number k. 
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The dependence of the growth rate on the wave number 
for ^ = 0.1,0.5,0.8. 



We analyze the growth rate 7 of small perturbations as a function 
of wave number k. When k < kc = jJ. the growth rate is posi- 
tive and any small perturbations will grow with time. For wave 
number k > kc, the growth rate 7 is negative. Consequently, 
the small perturbations always decay. It can be shown that, for 
wave number k = kc the growth rate is zero 7 = 0. As the wave 
number approaches zero k — > or infinity k — > 00 the growth 
rate approaches 7 = 1 and small perturbations grow with time. 
As the wave number approaches unity k — 1 the growth rate is 
determined by the value of jj, 



7 ^ 



2n 



M + 1 

For /i = the growth rate is zero 7 = 0. For = 1 the growth 
rate is unity 7 = 1. Since the growth rate is zero at the critical 
wave number and approaches to unity in the limit of zero wave 
number, there must exist a wave number in between that maxi- 
mizes the growth rate. Figs. (1-4) show the growing rate, 7, as 
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a function of the wave number k, for various values of the vortex 
mass fi. As the value of /u increases, the corresponding growth 
rate increases. 

Conclusion 

Thus, in the present work we show that under some condi- 
tions flux avalanche may occur in superconductor sample, which 
takes into account the inertial properties of the vortices. It has 
been noticed that a detailed theoretical study of this problem 
will be presented in our further papers. 
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Fig. 2-4. The dependence of the growth rate on the wave number for /i = 0.1, 0.5, 0.8. 
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